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^ : ABSTRACT 

We study the solution of viscous accretion disks around rotating com- 

, - pact/central object having hard surface i.e. neutron star, strange star and 

^ ' any other highly gravitating objects. We choose pseudo-Newtonian approach 



to describe the relativistic accretion disk. For this purpose, a new pseudo- 
Newtonian potential is established which is applicable to describe the rela- 
tivistic properties of star and its disk. As we know, the Hartle-Thorne met- 
ric can describe geometry of star as well as the space-time out-side of it, 
we use this metric to establish our potential. Our potential reproduces the 
marginally stable orbit exactly as that of general relativity. It also reproduces 
the marginally bound orbit and specific mechanical energy at the marginally 
stable orbit with at most 4% and 10% error respectively. Using this poten- 
tial we study the global parameter space of the accretion disk. Thus, we find 
out the physical parameter regime, for which the stable accretion disk can be 
formed around gravitating object with hard surface. We also study, how the 
fluid properties get changed with different rotations of the central star. We 
show that with the change of rotation to the central object, the valid disk 
parameter region dramatically changes. We also show the effect of viscosity to 
the fluid properties of the disk. Subsequently, we give a theoretical prediction 
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of kHz QPO, at least for one out of a pair, for a fast rotating compact object 
as 4U 1636-53. 

Key words: accretion, accretion disk — stars: rotation — hydrodynamics 
— gravitation — shock waves — X-rays: binaries 

1 INTRODUCTION 

For over last three decades, fluid dynamical study of accretion disks around compact object 
is indeed a topic extensively discussed, mostly for black holes (Abramowicz et al. 1988, 
1996; Abramowicz & Kato 1989; Chakrabarti 1990, 1996a,b,c; Narayan & Yi 1994; Narayan 
et al. 1997, 1998; Kato et al. 1998; Sponholz & Molteni 1994; Peitz & Appl 1997; Gammie & 
Popham 1998; Popham k Gammie 1998; Miwa et al. 1998; Lu k Yuan 1998; Manmoto 2000) 
and a few for neutron stars (Chakrabarti & Sahu 1997; Mukhopadhyay 2000, 2002a; Popham 
k Sunyaev 2001). Though in some of the cases authors use full general relativistic approach, 
in most of the cases the approaches are pseudo-Newtonian. In the pseudo-Newtonian method, 
they use simple non-relativistic equations to describe the disk with the appropriate choice 
of modified gravitational force which takes care about the geometry of back-ground space- 
time. There are a number of basic differences between the accretion disk around black hole 
and any other gravitating objects namely neutron star, white dwarf, strange star and other 
gravitating stars particularly at the inner edge of the disk. (1) While the black hole has event 
horizon, other gravitating objects mentioned above have hard surface. Thus, whereas in the 
case of other gravitating object close to the surface the matter speed must be subsonic, for 
black holes the matter speed is supersonic. (2) In case of other gravitating object that has 
hard surface, the matter comes into contact with its outer surface. Thus, there is an obvious 
effect of dragging to the exterior space-time due to its rotation. In terms of general relativity, 
corresponding space-time metric should be the solution for interior as well as exterior to 
it. Therefore, it is necessary to consider the 'dragging effect' of rotating gravitating star 
(particularly for neutron star, where the relativistic effect is very important) to the space- 
time. In case of black hole, which does not have hard surface, this effect comes to the 
space-time automatically through the 'Lense-Thirring drag' from the Kerr geometry which 
does not describe any physics interior to the horizon. (3) Further, the overall temperature 
as well as the luminosity of the disk is expected to be high in the cases of other gravitating 
star particularly for neutron star with respect to that of black hole (Mukhopadhyay 2000, 
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2002a) for a particular accretion rate. If the same parameters are chosen for the infalhng 
matter towards black hole and other gravitating objects separately, at same accretion rate, 
the maximum temperature in the disk around black hole will be lower at least a few factors 
times. Thus, in the disk around neutron star, expected generation/absorption of nuclear 
energy is higher than that of black hole (Mukhopadhyay 2002a) . As there is a clear distinction 
between accretion disks around black hole and other gravitating objects, hereinafter, we will 
call the class of gravitating objects having hard surface as GOHS to distinguish from black 
hole. GOHS can be meant as neutron star, strange star, white dwarf and any other highly 
gravitating objects where disk may form. As the radius of white dwarf is very large, the 
general relativity is not much important in the accretion disk around it. But our discussions 
will be generally true for disk around any kind of GOHS. 

Mukhopadhyay (2002a) has studied the accretion disk around weakly magnetized slowly 
rotating neutron stars, and found that the temperature of the inner disk is significantly 
higher than the case of black holes. Prasanna & Mukhopadhyay (2003) have described the 
accretion disk with the inclusion of 'Coriolis' effect. They have shown that the rotational 
effect of compact object to the accretion disk can be approximately reproduced with the 
use of Coriolis effect. They have done perturbative analysis and shown that the self-similar 
solutions are well behaved for most of the parameter space. As no cosmic object is static, 
in this paper, our interest is to include the rotation of the neutron star, globally speaking 
rotation of GOHS, in a more proper manner and see its effect onto the disk. The recent 
observation tells about the evidence of millisecond pulser (Strohmayer & Markwardt 2002) 
of frequency 582Hz for 4U 1636-53 which is the accreting LMXB. Therefore, it is important 
to consider the rotational effect of the GOHS to the disk. We approach our study in pseudo- 
Newtonian manner. So far, there is no pseudo-potential in the literature which can describe 
the accretion disk around GOHS. Thus, first of all we should establish such a potential useful 
for our work. 

Shakura & Sunyaev (1973) first initiated the accretion disk modeling using simple New- 
tonian potential meant for non-rotating black holes. As the relativistic effects are extremely 
important near the compact objects, this potential can not describe the essential inner prop- 
erties of the disk. Paczyhski & Wiita (1980) modified this Newtonian potential in conformity 
with the Schwarzschild geometry, which can naturally reveal approximately the properties 
of disk around non-rotating black holes and GOHSs, even that of the inner most part of the 
disk. The potential can reproduce the radius of marginally stable orbit (xg) and marginally 
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bound orbit (xj,), mechanical energy per unit mass at the last stable circular orbit (Eg) and 
viscous efficiency of the disk. The last two parameters agree within a 10% error but the first 
two cases exactly match with that of the Schwarzschild metric. Nowak & Wagoner (1991) 
proposed another potential for accretion disk around non-rotating black holes. This potential 
can also mimic approximately most of the properties of the disk governed by Schwarzschild 
geometry. However, angular and epicyclic frequency of the disk can be best analyzed with 
this potential. Artemova et al. (1996) proposed a couple of pseudo-potentials to describe 
accretion disk for rotating and non-rotating black holes. These potentials are well analyzed 
by a number of astrophysicists later (e.g. Lovas 1998; Semerak & Karas 1999). Very recently, 
Mukhopadhyay (2002b) and Mukhopadhyay & Misra (2003) prescribed some new potentials 
with which properties of accretion disk can be very well described in Kerr geometry. The 
more interesting fact lies in the methodology adopted by Mukhopadhyay (2002b) to describe 
a potential in the accretion disk, which can be used to derive the pseudo-potential for any 
metric according to the physics concerned. The potentials proposed by Mukhopadhyay & 
Misra (2003) can be used for the time-dependent simulation of accretion disk. But all the 
potentials are described either for the fluid dynamical study of disk around black hole or the 
study of temporal effect around compact objects. As the fluid properties of inner accretion 
disk are very much influenced by the nature of central gravitating object as mentioned above, 
it is essential to establish a pseudo-potential which can describe an accretion disk around 
GOHS. Although in case of black hole, the space-time can be described by Kerr geometry, 
but in the case of GOHS with slow rotation, the corresponding metric which can continu- 
ously describe the space-time inside and outside the star is Hartle-Thorne (hereinafter HT, 
Hartle & Thorne 1968). Thus, we should estabhsh our potential according to HT metric. 

Thus, we shall consider the general solutions for fluid disks around rotating GOHS. In 
order to consider a solution procedure, we will follow Chakrabarti (1996b) where the global 
solution for transonic accretion flows is studied around rotating black holes but in the case 
of weak viscosity. Here, we will consider the set of disk equations in terms of general viscous 
prescription and the accretion around GOHS. As we want to concentrate upon the features 
of the inner region, we shall restrict the study to 'sub-Keplerian' flow. We study, how does 
the rotation of GOHS affect on the global parameter space of accretion disk. We know that 
for the case of non-rotating compact object, the stable accretion disk forms only for a certain 
set of physical parameter. Wc would like to check how this parameter region gets affected, 
how does the steady disk structure change, with different values of the rotation to GOHS. 
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Once we have the clear picture about the parameter space, we will concentrate upon the 
fluid properties of the disk for different viscosity and how does the viscous fluid properties 
change with different rotations of GOHS. 

Therefore, we are going to present a complete set of work on accretion disk around GOHS. 
In order to do so, in the next section, we will derive the required pseudo-Newtonian potential 
and establish its efficiency. In §3, we will present the basic equations of viscous accretion 
disk around GOHS. In §4 and §5, we will analyse the parameter space of the accretion disk 
around GOHS and corresponding properties of accreting fluid respectively. Finally we shall 
present our discussion and overall summary in §6. 



2 PSEUDO-NEWTONIAN POTENTIAL AND ITS EFFICIENCY 

Here, following Mukhopadhyay (2002b), we will derive the pseudo-Newtonian potential for 
accretion disk modeling around GOHS according to HT metric. The Lagrangian density for 
a particle in the HT geometry (neglecting the higher quadruple terms) at the equatorial 
plane (^ = 7r/2) can be written as 

/ 2M 2J2\ .„ / 2M 2J^Y^ 2 2 -,2 4J,. 
2^ = -(l- — -— — + — j r^ + rV--0t, (1) 

where over-dots denote the derivative with respect to the proper-time r and J denotes the 
angular momentum of GOHS. As the metric is valid only for slowly rotating stars, the 
rotation parameter J is restricted as J < 0.5. The geodesic equations of motion are 

2M 2J2\ . 2J ■ 

E = constant =1 t H 0, (2) 

2J ■ 

A = constant = r^</) i. (3) 

r 

For the particle with non-zero rest mass, g^uP^P'^ = —fn? (where is the momentum of the 
particles and g^j_i, is the metric). Replacing the solution for t and from and into 
gives a differential equation for r 



dr J \ g(r)r^ g(r)r^^ r^ J V5'(^) 



1 (j((^rjr'^ g{r)r' I 
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where, g{r) = (l - ^ + h{r) = (l - ^ - Here, ^ can be identified as an 

effective potential for the radial geodesic motion. The conditions for circular orbits are 

*=0. - = 0. (6) 

Solving for E and A from (0) we get 

E = X{r), (6) 



X(r) 



\ [-2J3m2 + Jmh-'^{-Am + 3r) + ^m^ibJ"^ + Mr^){2.P + r^{-2M + r))^ 



(7) 



where, 

X(r) 



2rm' + Jm\AM - 3r)r'^ 



-^m4(5J2 + Mr3)(2J2 + r3(-2M + r))2 
{(-12jWr2 + 2JW(9M - 7r)r^ + m^M{3M - r)r^ 



+ 6Jrym4(5J2 + Mr^){2.P + r^{-2M + r))2)/ 
(36J^ - r*^(-3M + r)2 + 12jV2(-3M + 2r))}^/l 



(8) 



Now as standard practice, we can define the Keplerian angular momentum distribution. 



A;^ = -|. Therefore, corresponding centrifugal force in HT geometry can be written as 

^4^3/ 4J2 M^.3\2 



X2 



Fr. (9) 



[-2Pm? + Jm?r^{-Am + 3r) + ^m^{5P + Mr^){2P + r3(-2M + r))2]2 

Thus from above, can be identified as the gravitational force of rotating GOHS at the 
Keplerian orbit in an equatorial plane. If we choose x = r/M [as G = c = 1) and m = 1, 
above expression reduces to 

x^(-4j2 + x^)^ 



[-2J3 + Jx3(-4 + 3x) + ^(572 + x3)(2J2 + x3(-2 + x)) 



212 



(10) 



Also the above expression reduces to Paczyhski-Wiita (1980) form for J = 0. Thus we 
propose that (fTUj) is the most general form of the gravitational force corresponding to the 
pseudo-potential in accretion disk around GOHS at equatorial plane. 
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2.1 Comparison of the Results for Hartle-Thorne geometry and 
Pseudo-Potential 

Now we will establish the validity of this potential. We would like to check, whether this 
potential can reproduce the values of marginally bound {xb), marginally stable (xg) orbits 
and mechanical energy per unit mass at Xs (Es), as same as HT geometry or not. 

In terms of the peudo-potential, one can write the radial velocity (f), angular velocity 
and angular momentum per unit mass (A) for a Keplerian orbit in an accretion disk as 
(Lynden-Bell 1969, Pringle & Rees 1973, Shakura & Sunyaev 1973, Novikov & Thorne 1973, 
Paczyhski & Wiita 1980), 



v = Jx^, ^^ = J--, X = Jx^-. (11) 
V (ix ' \ X dx'' \ dx 

At the marginally bound orbit, mechanical energy E reduces to zero and we get 

^ + y = :^^ + I/ = 0, (12) 
2 2 dx ^ ' 

where, V = J Fdx, which can be used to calculate Xb for the potential. For the last stable 
orbit (xs), dX/dx = 0, for our pseudo-potential which gives 



x^ 



-8J^ + 2J2(5 _ 2x)x^ + {x- 2)x^ -2407** + (x - 6)x^ + 8JV(12 + 5x) 



9 I Q 74^3/ 



0. (13) 



+2J^x%20x - 33) + (24 J3 + 127x^)^(572 + x3)(2J2 + - 2)x3)2 

The solution of p3|l gives the location of the last stable circular orbit (xs). For any J, the 
value of Xg computed from the above equation (fT^. matches exactly with the radius of last 
stable circular orbit in HT geometry. In Table-1, 2 and 3, we list Xb, Xg and Eg respectively 
for the potential V and HT geometry for various values of J. 
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Table- 1 
Values of Xb 



J 





0.1 


0.2 


0.3 


0.4 


0.5 


V 


4 


3.7862 


3.5619 


3.3255 


3.0762 


2.8135 


HT 


4 


3.7954 


3.5810 


3.3557 


3.118 


2.8689 


J 





-0.1 


-0.2 


-0.3 


-0.4 


-0.5 


V 


4 


4.2042 


4.4 


4.5885 


4.77 


4.9462 


HT 


4 


4.1957 


4.3838 


4.5684 


4.7397 


4.9089 



Table-2 

Valuers of ,Ts 





J 


0.1 


0.2 


0.3 


0.4 0.5 




y & HT 


6.0 5.6648 


5.3107 


4.9339 


4.5299 4.0934 




J 


-0.1 


-0.2 


-0.3 


-0.4 


-0.5 




V km 


6.0 6.3189 


6.6240 


6.9170 


7.1994 7.4724 


Table-3 
Values of Eg 


J 





0.1 0.2 


0.3 


0.4 


0.5 


V 


-0.0625 


-0.0664 - 


-0.0711 


-0.0769 


-0.0843 


-0.0943 


HT 


-0.0572 


-0.0607 - 


-0.0649 


-0.0701 


-0.0767 


-0.0857 


J 





-0.1 


-0.2 


-0.3 


-0.4 


-0.5 


V 


-0.0625 


-0.0592 - 


-0.0563 


-0.0538 


-0.0516 


-0.0497 


HT 


-0.0572 


-0.0542 - 


-0.0517 


-0.0494 


-0.0474 


-0.0457 



From Tabic- 1, it is clear that for all values of J, V can reproduce the value of Xb in very 
good agreement with general relativistic (HT) results. The maximum error in xi, is ~ 4%. 
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Table-2 indicates that V reproduces Xg exactly as that of HT geometry. Table-3 shows a 
maximum possible error in Eg is ~ 10%. Thus, the potential V will produce a slightly larger 
luminosity than the general relativistic one in the accretion disk for a particular value of 
J. It is to be noticed that for counter-rotating GOHS, that is for retrograde orbits, the 
errors are less than those of co-rotating ones. Eventually, we can tell that our potential can 
describe approximately all the phenomena of that of HT geometry and can claim it as a 
good pseudo-potential to describe the relativistic accretion disk around the rotating GOHS, 
particularly close to the equatorial plane. Thus, for a thin disk or vertically averaged disk, 
this is an ideal pseudo-potential. 



3 BASIC EQUATIONS OF ACCRETION DISK AND THEIR 
CONSEQUENCES 

Here, throughout our calculations, we express the radial coordinate in unit of GM/c^, where 
M is mass of the compact object, G is the gravitational constant and c is the speed of light. 
We also express the velocity in unit of speed of light and the angular momentum in unit of 
GM/ c. The equations to be solved are given below. 

(1) Equation of continuity: 

- 47rxEv = M, (14) 

where E is the vertically integrated density and M is the accretion rate. Following Mat- 
sumoto et al. (1984), we can calculate the E as 

S = InPeh{x), (15) 

where Pe is the density at equatorial plane, h{x) is the half-thickness of the disk and /„ = 
(2n+i)! ' "^li^re n is the polytropic index. From the vertical equilibrium assumption, the half- 
thickness can be written as 

h{x) = CsX^/^F-^/^ (16) 

where Cg is the speed of sound. 

(2) Radial momentum balance equation: 

dx p dx x^ ' 
where the flow is chosen adiabatic with the equation of state to be P = kp^ and = 

(3) Azimuthal momentum balance equation: Here, we will follow Chakrabarti (1996a) to 
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express the viscous dissipation in terms of shear stress Wxip- Because of this viscous 
dissipation, the angular momentum varies in accretion disk. Thus, Wxif, = —a{In+iP + 



Inv'^p)h{x) and Q~ 



x<p 



7] is coefficient of viscosity and a is Shakura-Sunyaev (1973) 



viscosity parameter. Thus the equation to be 



1 



d 



dx ph{x)x dx 



x^a 



(18) 



(4) Entropy equation: According to the mixed shear stress (Chakrabarti 1996a), = 
—a{In+iP + Inv"^ p)h{x)x^. For simphcity, we also consider the heat lost is proportional to 
the heat gained by the flow. Thus the equation to be 



Y^vT— = '"^^^^ (—- Ti-^\ =Q+ -Q- 
dx Ts — 1 \dx p dx J 



(19) 



where s is entropy density and / is cooling factor which is close to and 1 for the flow with 
efficient and inefficient cooling respectively. Following Cox & Giuli (1968), we can define 
Ti - /? ^ ^ , (4 - 3/3)2(7 - 1) ^ pkT/pmp 



T 



1 + 



Fi =/3 + 



(20) 



4-3/5' ' /3 + 12(7-1) (1-/5)' aTy3 + pkT/pmp' 

Here, P (ratio of gas pressure to total pressure) close to for radiation dominated flow 
(highly relativistic flow) and close to 1 for gas dominated flow. Now combining ()14|) . ()17p . 
(UHl) and (Uni), we get 



dv_ _ f\{x,v,c,) 

dx f2{v,Cs) 

where, 

fl{x,V,Cs) = 



(21) 



\2x 2F dx) ^ ' 
2ad 



x-' 



f2iv,Cs) 



1 



(37 — l)xv 
2c? 



+ 



av 



X 



2A 



A + 



3_ 

2x 



1 dF" 
2F~dx . 



4aA 



t;(37- 1) 
vcliV, - 1), 



aA-c^(Fi - 1) + 7t; 



1 

V 



4aA 



(37 - l)f 



(22) 
(23) 



(37- l)f2_ 

and A = 7(F3-l)/a(^f + 

As far away from the black hole, v < Cs and close to it, v > Cg, there must be an in- 
termediate location, where the denominator of (j?H) must be vanished. To have a smooth 
solution, at that location, numerator has to be zero. This location is called the sonic point 
or critical point (xc)- The existence of sonic location plays an important role in accretion 
phenomena. Although, in case of accretion flow of matter around black hole sonic point must 
exist, for that around GOHS (say, neutron star) it is not always necessary (Mukhopadhyay 
2002a). When shock forms in an accretion disk around a neutron star (Chakrabarti & Sahu 
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1997, Mukhopadhyay 2002a), sonic points play an important role. From the global study 
of the sonic points in an accretion disk, one can understand about the stability of physical 
parameter region, which we will discuss in the next section. 

Now, when the sonic point exists in an accretion disk around GOHS, f2{vc,Csc) = at 
X = Xr- Thus at that location Mach number can be written as 



. ^±«^. (24) 

where 

A = „V(r3-i) + 7(r. + i)-^^^f^. 

37-1 V /„ 

37-1 /n ^ ^ 

Subsequently, at a; = Xc, fi{xc,Vc, Csc) = 0. Thus, using we can eliminate Vc from 
fi{xc, Vc, Csc) = and get an algebraic equation for Cgc and Xc, which we can solve to find out 
sound speed at sonic location. To find out ^|c, we apply I'Hospital's rule to (j2H). In order to 
understand the fiuid properties in accretion disk, we have to solve (PT|) with an appropriate 
boundary condition. Also, integrating p8|) . we get the angular momentum of the accreting 
matter as 

where M is the Mach number of the fiow and Aj„ is the angular momentum at the inner 
edge of accretion disk. 

Now integrating (fT^ and p7|) . we can write the entropy and energy of the fiow at sonic 
point as 

7-1 

(27) 




,,3 -^c 

(7-1) '~ 



(28) 



where K = {I Fdx)\e and q = It can be mentioned that, disk entropy and accretion 

rate are related by a simple relation as = {''yK)^M. While M is a conserved quantity 
for the particular accretion fiow, is not, as it contains which carries the entropy 
information that is not conserved in a dissipative system. As a boundary condition for the 



12 Mukhopadhyay & Ghosh 

particular flow, we have to supply the sonic energy E^. Then from (j^Hj) . we can find out the 
sonic location Xc- Therefore, knowing Xc one can easily find out the fluid velocity and sound 
speed at the sonic point from /i = /2 = with the help of for a particular accretion 
flow. These have to be supphed as further boundary conditions of the flow. 

Now let us come to the issue of the formation of shock. Mukhopadhyay (2002a) showed 
that in accretion disk around neutron star, double shock is very natural in certain physical 
situations. Here, whenever we mention the shock, we will mean Rankine-Hugoniot shock 
(Landau & Lifshitz 1987). If we generalise the conditions to form a shock in an accretion 
disk given by Chakrabarti (1989) for rotating GOHS, we get 

^Mlc^^ + ncl^ = + ncl_, (29) 

( 27 , . ,,2\ ( 27 ^ 



M+>M^, (31) 
where 

^3/2 

7W = Mc^("+^)-^. (32) 



Here, subscript '— ' and '+' indicate the quantities just before and after the shock respec- 
tively, Xs indicates the shock location and v = ^E^- Fi'O'^ (j29j) - (j32|) . it is very clear that 
except the entropy expression, all remain unchanged with respect to the case of a non- 
rotating central object. Also from ()30|) and ()32j). we get the shock invariant quantity as 

^ _ (ivfc + (37 - m^y _ {it + (37 - m~y 

M2(7-l) + 2 M2(7-l) + 2 ' ^ ^ 

which remains unchanged with respect to a non-rotating case. If only all the conditions, 
(f^ - (p?T|l and pHjl are simultaneously satisfied by the matter, shock will be generated in an 
accretion disk. 



4 GLOBAL ANALYSIS OF THE PARAMETER SPACE 

One of our aim is to check, how does the rotation of GOHS affect the disk parameter 
region known for non-rotating case. More precisely, we will check, how does the rotation 
of GOHS affect the sonic location, structure as well as the stability of disk. As mentioned 
in §1, we will follow the methodology adopted by Chakrabarti (1996b) for global solution 
and Mukhopadhyay & Chakrabarti (2001) for stability analysis. In Fig. 1, we show the 



Global solution of accretion disk around rotating compact objects 13 
variation of disk entropy as a function of sonic location. The intersection to the horizontal 
line (which indicates a constant entropy curve) with all the curves indicate the sonic points 
of the accretion disk for that particular entropy and rotation to the GOHS. It is clearly seen 
that at a particular entropy, if the co-rotation of GOHS increases, the inner sonic points shift 
to more inner edge of the disk, which are comparatively unstable region (because, beyond of 
Xg, disk orbit is not stable and decrement of the value of Xg is slower than that of Xc with J). 
For J = 0.5, the inner sonic point disappears, as we see that dashed curve terminates on the 
radial axis at ~ 6 and it does not turn up to give rise the inner sonic point. The reason 
behind it is that, for higher the co-rotation, the angular momentum of the system will become 
higher and radial matter speed will be unable to overcome the centrifugal barrier close to 
the GOHS and that will make the disk unstable at the inner edge. For the counter-rotating 
cases, if J increases in magnitude for a particular entropy, disk will become more stable, as 
the radial speed will be able to overcome the centrifugal barrier at a larger radius compared 
to a co-rotating case and matter may attain supersonic speed at a larger distance, which is 
more stable region of the disk. However, for a high counter rotation, angular momentum of 
the system will become too low and the accretion will tend to 'Bondi-hke flow' (Bondi 1952) 
with the existence of single sonic point. 

It is seen that for a = — 0.1, the variation of curves are almost same and only significant 
changes come into the picture for a > 0.1. Also for the lower viscosity, if the co-rotation is 
high (J is high), the possibility to have the all three sonic points in the disk is high. As the 
viscosity increases, the possibility for the existence of all three sonic points diminishes. In 
order to preserve the existence of sonic point for the high viscous cases, J must reduce to 
zero and subsequently, more negative. As a whole, from Fig. la-c we can conclude that for 
lower a, multiple sonic points may exist, but for high a, there is only one sonic point which 
is inner one in the accretion disk. Similar features come out from Fig. 2 but in terms of 
energy of the sonic point in accretion disk. As we know that the observed angular frequency 
for the candidate 4U 1636-53 is 582Hz, which is equivalent to J = 0.2877, one of the value 
of J chosen for our discussion is 0.2877. For both Fig. 1 and 2, sonic points with negative 
slope of the curve indicate the locations of 'saddle-type' or 'X-type' sonic point and positive 
slopes indicate the 'centre-type' or '0-type' sonic point. Thus, the rotation of GOHS and 
viscosity of the disk are the important factors to the formation and location of sonic points, 
which are related to the structure of accretion disk. 

Figure 3 indicates the variation of energy with entropy at sonic locations, when J is 




Figure 1. Variation of sonic entropy as a function of sonic location for various values of J when (a) a = 0,/ = 0, (b) 
a = 0.4,/ = 0.5, (c) a = 0.8,/ = 0.5. Central solid curve indicates non-rotating case (J = 0), while the dotted and dashed 
curves are for J = 0.2877, 0.5 respectively and long-dashed and dot-dashed curves are for J = —0.2877, —0.5 respectively. The 
horizontal line indicates the curve of constant entropy of 7 X 10~®. The other parameters are A = 3.3, 7 = 4/3. 



parameter for different viscosity of tfie fluid. It is seen that, witfi the increase of J (co- 
rotation), inner sonic points of the disk shift towards the lower entropy region keeping the 
outer sonic points unchanged and the disk becomes unstable. We know that the shock in 



accretion disk can be formed only if there is a possibility of matter in the outer sonic point 
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Figure 2. Same as0 but sonic energy is plotted in place of entropy. The liorizontal line indicates the curve of constant energy 
of 0.005. 

branch of lower entropy to jump to the inner sonic point branch of higher entropy. As 
because, with the increase of J, inner sonic point branch itself shifts to the lower entropy 
region which is comparatively unstable, the shock as well as disk also become unstable. For 
J = 0.5, the inner sonic point branch disappears totally and only outer and 0-type sonic 
point branches remain. Clearly, for a highly rotating GOHS, the possibility of formation of 
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shock reduces as well as the shock becomes unstable, as the inner sonic point branch itself 
disappears or tends to disappear. The physical reason behind it is that, with the increase of 
value of J, the angular momentum of the system increases, which helps the disk to maintain 
a high azimuthal speed of matter upto a very close radius to the GOHS, as a result the radial 
speed of the matter can overcome the corresponding sound speed at a very inner radius only. 
On the otherhand, the inner edge of the accretion disk is comparatively less stable, as the 
entropy decreases at lower radii. As the sonic points form at a more inner radii, with the 
increase of J, disk tends to an unstable situation. 

In case of the counter-rotation of GOHS, as J increases in magnitude, the inner sonic 
point branch shifts towards a higher entropy region, and the disk as well as the shock 
become more stable. The physical reason to have a more stable inner sonic point branch 
for the counter-rotating case is the following. As the counter rotation of central object 
increases, that is the angular momentum of the system reduces, the matter falls freely 
towards GOHS at larger radii and it becomes supersonic at a comparatively outer edge of 
the disk. Thus, the disk as well as shock become more stable. However, for a very high 
counter-rotation of GOHS, this inner sonic point branch itself tends to merge to the outer 
sonic point branch, and the entire system shift towards the outer side which is more stable. 
As in case of higher counter-rotation, the angular momentum of the system becomes very 
small, the corresponding centrifugal pressure onto the matter becomes insignificant, as a 
result the possibility of shock diminishes again. Thus, if the angular momentum of GOHS 
increases or decreases significantly, the shock wave in accretion disk becomes unstable and 
the disk itself tends to an unstable situation. In Fig. 3, 'C' is the point, where the flow can 
pass through the inner and outer sonic point simultaneously. The parameters around C are 
important to study the shock in accretion disk. 

As the viscosity of flow increases, the possibility of intersection between inner and outer 
sonic point branch decreases. With the increase of viscosity, more regions containing the X- 
type sonic points become 0-type. The outer X-type sonic points recedes in further out and 
inner one proceeds in further inwards. Also for higher viscosity, the outer sonic points may 
no longer remain as X-type and only the inner sonic points may exist. Thus, the possibility 
of shock will reduce. 
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5 FLUID PROPERTIES OF ACCRETION DISK 

Now we will discuss the behaviour of viscous accreting fluid around rotating GOHS. In 
the previous section, we have already conflrmed that the rotation significantly affects the 
global parameter space of accretion disk. Here, our intention is to see, how the various fluid 
dynamical results of an accretion disk are affected for various values of J and a. For a certain 
choice of physical parameter set, if the existence of sonic radii are possible in the accretion 
disk around GOHS, the flow structures can be studied for the corresponding choice of sonic 
radii (or energy) of the accretion flow. 

In Fig. 4a, we show the variation of Mach number for three different values of J at a 
particular viscosity. As the viscosity is low, the rate of energy momentum transfer in matter 
is less. It is known that for an accretion around neutron star, in a certain cases, matter may 
be always subsonic (Chakrabarti & Sahu 1997; Mukhopadhyay 2002a) and at the stellar 
surface matter speed reduces to zero. Owing to this fact, here such a situation is considered 
where the speed of the accreting fluid is low, which accounts for the high residence time 
of matter in the disk. The angular momentum of the disk is high at x ~ 10 (see Fig. 4d) 
and thus radial matter speed is low. Subsequently, though the disk angular momentum goes 
down, as the matter comes closer to the stellar surface, it starts to decelerate and the radial 
matter speed still remains low. This results to a higher the possibility of cooling, and thus 
cooling factor / is chosen as 0.1. As the energy momentum transfer rate of matter is low, the 
accretion rate is chosen intermediate. We have chosen a standard mass of GOHS as 2Mq. 
If the GOHS is chosen to be non-rotating, angular momentum of the system is less and the 
centrifugal barrier is minimum. With the increase of J, this centrifugal barrier increases, and 
it becomes very high particularly for J = 0.5. Mach number proflles clearly indicate that 
at the surface of GOHS, accreting fluid stops. Also from Fig. 4b and 4c, it comes out that 
sudden deceleration of the accreting fluid gives rise to a sudden enhancement of density and 
temperature close to GOHS. Here, the temperature is thought to be cooled down by a factor 
of 1/30 through the inverse-Compton effect in the radiation pressure dominated relativistic 
flow. Thus, we choose P ~ 0.03. Nevertheless, the temperature is still very high, which 
was also reported earlier (Mukhopadhyay 2002a) in the context of non-rotating neutron 
star. The density proflles are plotted in unit of qT]^- Figure 4d shows the ratio of sub- 
Keplerian to Kcplerian angular momentum variation of the accreting fluid as a function of 
radial coordinate. It is seen that, at the boundary between Keplerian and sub-Keplerian 
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disk, Afc/A 1, and from that radius (xk) we start our calculation as our discussion is 
for sub-Keplerian accretion disk. From the A^/A profile, it is again clear that for higher J, 
centrifugal barrier becomes stronger. Far away from GOHS, the effect of centrifugal barrier 
is weak, but as the matter approaches towards GOHS, it starts to dominate and at x ~ 10 
it becomes stronger. Then again it goes off. However, at very close to the GOHS, again this 
centrifugal effect comes into the picture as the matter falls into the strong rotational field of 
GOHS. This effect appears very prominently for high value of J, which is here chosen to be 
0.5. Actually, this effect comes into the picture, once the matter enters into the corresponding 
ergosphere of HT geometry. But the important point to be noted that if the outer radius of 
GOHS is greater than the radius of ergosphere, this feature is not possible to occur. 

In Fig. 5, we show the results of accretion phenomena for 4U 1636-53, whose angular 
momentum J = 0.2877. The examples are shown, where the sonic point (s) exist in the flow 
and shock forms in the accretion disk around 4U 1636-53 for different viscosity parameters. 
Figure 5a shows the variation of Mach numbers as a function of radial coordinate. It reflects 
that, for a small a (~ 0) only shock is possible at x = 19.02 in the accretion disk. As a 
increases, the energy momentum transfer rate increases, rate of inf ailing the matter enhances 
and velocity of the accreting fluid becomes high at the inner edge of the disk. But, because 
of the inner hard surface, matter has to stop at close to GOHS and naturally another shock 
forms at the inner edge of the accretion disk. The shock locations for a = 0.02 and 0.05, 
are at x = 14.27,4.38 and x = 15.37,4.31 respectively. If the viscosity is high, the cooling 
factor is intermediate (see Fig. 5 caption). Also, the angular momentum of matter at the 
sonic point is chosen as, Ac = 3, and for the radiation dominated, inverse- Comptonised flow, 
P is chosen to be 0.03. Figures 5b and 5c show the variations of density and temperature. 
The regions, where velocity decelerates abruptly, temperature and density also jumps up in 
a significant manner. Here also, the inverse- Comptonised temperature is very high as the 
cases of Fig. 4. 

6 SUMMARY AND DISCUSSION 

Here, we have studied the viscous accretion phenomena around rotating gravitational object 
with hard surface (GOHS), i.e., mostly around compact object like neutron star and strange 
star. We started with a proposal of a new pseudo-Newtonian potential by means of Hartle- 
Thorne geometry, as this metric can describe the space-time of interior as well as exterior 
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to a slowly rotating star. Wc have shown that this potential can describe all the essential 
relativistic properties like, radius of marginally boTind and stable orbit, specific mechanical 
energy etc., of accretion disk within 10% error. Then, along with this potential, we have 
described the accretion disk by using non-relativistic equations. Therefore, we can say, our 
method is a 'pseudo-general- relativistic'. Once the pseudo-Newtonian potential is proposed, 
we have applied it to the study of relativistic fluid properties of accretion disk around rotating 
GOHS. We have prescribed the set of basic equations for vertically averaged, thin, viscous 
accretion disk. Then, we have analysed the viscous parameter space globally. Subsequently, 
we have discussed about the properties of viscous fluid and the effect of rotation on the 
fluid behaviour. Examples are shown for the accretion disk around an observed candidate 
4U 1636-53, which is one of the fast rotating compact object with angular frequency 582Hz 
(equivalent J — 0.2877). As starting from a new pseudo-potential, we have discussed upto 
the solution of viscous accretion flow, this is one of the most self consistent paper on accretion 
disk. 

We have studied the stabihty of accretion disk. Most of the earlier analysis of structure 
and stability of the accretion disk have been done around non-rotating central object. As 
no cosmic object is static, it is very important to consider the effect of rotation, particularly 
for the discussion of an inner edge of the disk. We have found that, when the rotation to a 
GOHS is incorporated, the locations of sonic point (if any) get shifted to the unstable region 
and the shock (if any) gets affected in the disk, and thus the disk structure gets infiuenced. 
Therefore, we can conclude that the rotation should be incorporated in the study of an 
accretion disk. Any related conclusion is dependent on the rotation parameter of GOHS. 
The mentioned disk properties and phenomena not only get modified and/or shifted in 
location, sometimes disappeared completely. Prasanna & Mukhopadhyay (2003), worked on 
the stability analysis of accretion disk around rotating compact objects in a perturbative 
manner. They incorporated the rotational effect of central object in a disk indirectly by the 
inclusion of Coriolis acceleration term. Here, this rotational effect is brought directly from 
the metric itself. 

Prom the global analysis of sonic points, we have seen that for a higher co-rotation 
of GOHS, the disk becomes unstable for a particular angular momentum of the accreting 
matter. As the formation of shock needs a stable inner sonic point in the accretion disk, for 
the higher value of J, shock is unstable as the inner region of disk is unstable. On the other 
hand, for counter-rotating cases, the angular momentum of system reduces and the matter 
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falls more strongly to a GOHS. As there is no significant centrifugal barrier to slow down 
the matter in disk, the possibility of shock reduces again. Also the branch of inner sonic 
point merges or tends to merge to that of outer sonic point. Thus, we can conclude that the 
parameter region, where the shock is expected to form for non-rotating GOHS (Chakrabarti 
& Sahu 1997), is affected for rotating ones and the possibihty of shock is reduced for both 
the co-rotating and counter-rotating GOHS. Also for the higher viscosity, both the inner 
and outer sonic point branch merge each other and the possibility of shock reduces again. 
Actually, with the increase of viscosity, the region containing the X-type sonic points tends 
to become 0-type. The outer X-type sonic points recede in more outwards and inner ones 
proceed in inwards more. 

As the incoming matter slows down abruptly at the surface of GOHS as well as at two 
shock locations (if two shocks form), the overall density becomes higher for an accretion 
disk around GOHS compared to that for a black hole. Similarly, the temperature in the 
disk is higher for a GOHS. Thus, one can conclude that the accretion disk around a GOHS 
is very favourable for nucleosynthesis. We know that the accretion disk around a black 
hole is enough hot for nucleosynthesis (Mukhopadhyay 1999; Mukhopadhyay & Chakrabarti 
2000), which are different from that in star. As the density and temperature may be much 
higher around a GOHS, more efficient nucleosynthesis is expected particularly at an inner 
region of the disk. In early, we saw that the high temperature of accretion disk around a 
black hole is very favourable for the photo-dissociations and the proton capture reactions 
(Chakrabarti & Mukhopadhyay 1999). As the accretion disk around a GOHS is hotter than 
that around a black hole, even '^He which has a high binding energy, may dissociate into 
deuterium and then into proton and neutron. If we consider the accreting matter to come 
from the nearby 'Sun like' companion star, the initial abundance of ^ife in the accreting 
matter is supposed to be about 25%. Therefore, by the dissociation of this '^He, neutron 
may produce in a large scale, which could give rise to the neutron rich elements. Guessoum 
& Kazanas (1999) showed that the profuse neutron may be produced in the accretion disk 
and through the spallation reactions lithium may be produced in the atmosphere of the star. 
When the neutron comes out from the accretion disk by the formation of an outflow, in that 
comparatively cold environment, ^Li may be produced which can be detected on the stellar 
surface. In early, it was shown that the metalicity of the galaxy may be influenced when 
outflows form in the hot accretion disk around black holes (Mukhopadhyay & Chakrabarti 
2000). In case of the lighter galaxy, the average abundance of the isotopes of Ca, Cr may 
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significantly change. Also the abundance of lighter elements, like the isotopes of C, O, Ne, 
Si etc. may increase significantly. As the temperature of the accretion disk around GOHS is 
higher, the expected change of abundance of these elements and the corresponding influence 
on the metahcity of the galaxy is expected to be high. 

We know that, out of an observed pair of kilohertz QPO frequencies for a particular 
candidate, one is the oscillation due to its Keplerian motion (Osherovich & Titarchuk 1999). 
For 4U 1636-53, lower frequency is 950Hz (Wijnands et al. 1997) which may be the Keplerian 
one. According to our pseudo-potential, if wc calculate this Keplerian frequency for 4U 1636- 
53, it comes out to be 1027Hz for M — IAMq and xk — 18, which is not much over estimated 
with respect to the observed one. Thus, we can give a theoretical prediction of QPO, at least 
for one, out of a pair. Similarly, for other candidates, one can calculate Keplerian frequency 
according to our potential and compare with observation. The detailed theoretical study of 
QPO based on this present scheme will be pursued elsewhere. 
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Figure 3. Variation of sonic energy as a function of sonic entropy for various values of J when (a) q = 0, / = 0, (b) 
a — 0.4,/ = 0.5, (c) OL = 0.8,/ = 0.5. Solid curve indicates non-rotating case (J = 0), while the dotted and dashed curves are 
for .1 = 0.2877, 0.5 respectively and long-dashed and dot-dashed curves are for J = —0.2877, —0.5 respectively. O and I indicate 
the locus of outer and inner sonic point respectively, while C is the intersection of outer and inner sonic point locus. The other 
parameters are A = 3.3, 7 = 4/3. 
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Figure 4. Variation of (a) Mach number, (b) density, (c) temperature in unit of lO", (d) ratio of sub-Keplerian to Keplerian 
angular momentum of the accreting fluid as a function of radial coordinate. Solid, dotted and dashed curves are for J = 0, 0.1, 0.5 
respectively. Other parameters are a = 0.0001, / = 0.1, M = 2, M = 2, /3 = 0.03. 
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Figure 5. Variation of (a) Mach number, (b) density, (c) temperature in unit of lO" as a function of radial coordinate for the 
accretion disk around 4U 1636-53. Solid, dotted and dashed curves are respectively for (i) a = 0, / = 0, (ii) a = 0.02, / = 0.2, 
(iii) a = 0.05, / = 0.5. Other parameters for 4U 1636-53 are J = 0.2877, Ac = 3, M = 1.4, M = 1, ^ = 0.03. 



